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Based on the recently developed picture of an electronic ideal relativistic fluid at the Dirac 
point, we present an analytical model for the conductivity in graphene that is able to describe 
the linear dependence on the carrier density and the existence of a minimum conductivity. 
The model treats impurities as submerged rigid obstacles, forming a disordered medium 
through which graphene electrons flow, in close analogy with classical fluid dynamics. To 
describe the minimum conductivity, we take into account the additional carrier density 
induced by the impurities in the sample. The model, which predicts the conductivity as a 
function of the impurity fraction of the sample, is supported by extensive simulations for 
different values of S, the dimensionless strength of the electric field, and provides excellent 
agreement with experimental data. 
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Graphene has revealed a series of amazing properties, such as ultra-high electrical con- 
ductivity [l,2], ultra-low shear viscosity to entropy ratio {3], as well as exceptional structural 
strength, as combined with mechanical flexibility (4] and optical transparency {5]. Many of 
these fascinating properties are due to the fact that, consisting of literally one single carbon 
monolayer, graphene represents the first instance ever of a truly two-dimensional material 
(the "ultimate flatland" jsj). Moreover, due to the special symmetries of the honeycomb lat- 



fluid of massless chiral 

n 



tice, electrons in graphene are shown to behave like an effective Dirac 

quasi-particles, propagating at a Fermi speed of about c ~ 10 6 m/s |3|, |7j|. This configures 
graphene as an unique, slow-relativistic electronic fluid, where many unexpected quantum- 



electrodynamic phenomena can take place 



For instance, since electrons are about 300 



times slower than photons, their mutual interaction is proportionately enhanced, leading to 
an effective fine-structure constant a gr = e 2 jfivp ~ 1. As a result of such strong interac- 
tions, it has been recently proposed that this peculiar 2D graphene electron gas should be 
characterized by an exceptionally low viscosity/entropy ratio (near-perfect fluid), coming 
close to the famous AdS-CFT lower bound conjectured for quantum-chromodynamic fluids, 
such as quark-gluon plasmas js]. This spawns the exciting prospect of observing electronic 



pre-turbulence in graphene samples, as first pointed out in Ref. |3[ and confirmed by recent 
numerical simulations {9]. 

Some of the electrical properties of graphene are still not fully understood, such as the 
linear increase of the electrical conductivity with the number of charge carriers, the existence 
of a minimum conductivity (see Ref. 10J, and reference therein), and even the nature of 



the main scattering mechanism limiting the carrier mobility 



11] . In fact, classical transport 



theories, based on short-range scattering of electrons by impurities, predict that the electric 



conductivity in graphene should be independent of the carrier density 12[ . Recent works in 



vity m j 

HQ 



the field [131, M have shown that such linear dependence might be potentially explained by 
treating the impurities as screened Coulomb scatterers. Nevertheless, some measurements 
of the change in the electrical conductivity upon immersion of graphene samples in high- 
k dielectric media differ from this conclusion 
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16| . Here, we construct a model for 



describing the electrical conductivity in graphene by using a completely different approach, 
which is based on the recently developed picture of an electronic ideal relativistic fluid at 
the Dirac point. We demonstrate that, although this model is based on a semiclassical 
theory (it cannot take into account all quantum effects, e.g. Landau quantization, quantum 



hall effects, and quantum interference), it captures the main factors that contribute to such 
linear behavior and the appearance of a minimum conductivity. 

Since the most likely relevant limiting factor for the graphene conductivity is still subject 
of controversy, e.g. it can be due to random charged impurity centers 13[ or strong neutral 



defects that induce resonant scattering [15|, [16[, we will treat the impurities as hard-spheres, 
hindering the electron flow (scattering electrons), similarly to the way a disordered medium 
does in the context of fluid dynamics. The choice of hard-spheres is based on the experimen- 
tal results by Monteverde et al. ll|, which suggested that electrons seem to collide mostly 
with short range scatterers of the size of a few carbon-carbon interatomic distances, like 
voids, adatoms, etc. Since the relativistic fluid approach is derived from the quantum Boltz- 
mann equation (QBE) for graphene [l?}], a hydrodynamic description of the conductivity 
can be expected to apply as long as the QBE collision operator takes proper account of the 
Coulomb interactions between electrons. Therefore, once Coulomb interactions are included 
in the viscosity of the fluid, the conductivity (which in our case, unlike viscosity, is a prop- 
erty of the flow rather than of the fluid) becomes a function of the dissipation introduced in 
the system by the impurities, i.e. the electron-impurity interaction. 

Here, we treat graphene as a disordered medium and develop a model for its conduc- 
tivity, as a function of the impurity density describing the anomalous dependence of the 
conductivity on the carrier density and the minimum conductivity due to the carrier density 
induced by the impurities. The results are compared with experimental data yielding very 
satisfactory agreement. 



RESULTS 



Electronic Fluid in Graphene 



Our treatment is based on the hydrodynamic equations derived by Miiller et al 
based on the quantum Boltzmann equation for electrons in graphene. This analysis delivers 
the value of the transport coefficients, namely the fluid shear viscosity, which is an input 
parameter in our model. The hydrodynamic approach in graphene is valid when the inelastic- 
scattering rate due to electron-electron interactions dominates. This is the case at low 



doping, at high temperatures, and in moderate fields 



181 ] . However, to neglect electron- 



phonon interactions, we will have to stay at a moderately high temperature of around 



100K 



19]. In absence of magnetic fields, the quasiparticle distribution function, f s , evolves 



according to the quantum Boltzmann equation, 

df s 



at 



+ v s -Vf s + eE-V l :f s = -n[f s 



(1) 



where E is an external electric field, e the electric charge of the electron, Q[f s ] a collision 
operator that takes into account the electron-electron interactions, v s = s c k/\k\, k the 



wave vector, c the Fermi speed (~ 10 6 m/s), and the sign s, not to 



:>e confused with the 



entropy density, distinguishes between electrons (+) and holes (— 
the probability distribution function becomes the Fermi-Dirac distribution 

1 



201 ] . At equilibrium, 



f s {t, x, k) 



(2) 



e (sc\k\-/i)/k B T _|_ i 

where p is the chemical potential and T denotes the temperature. Thus, in the hydrodynamic 
limit, from Eqs. ([Q) and (T5]) one can derive the equations for the Dirac electron fluid in 
graphene: dp/dt + V ■ (pu) = 0, for charge conservation; de/dt + V ■ [(e + p)u] = 0, for 
energy density conservation and 



Pr 



du 
~dt 



XL UT) —* 

+ Vp + — -£ -7]V 2 u = pE 
c 2 at 



(3) 



for momentum conservation. Here, e is the energy density, p the pressure, p the charge 
density, u the drift velocity, p r = (e + p)/c 2 , and t] the shear viscosity. 

For the case of undoped graphene (p — 0), the presence of charge density is due to the 
thermal energy and can be described by, 

P = p th = e (^) . (4) 

However, when there are impurities, they can induce electric potentials on the graphene 
sample and a correction due to the chemical potential must be considered, 



P = pth§p{p/k B T) 



where $ p is a dimensionless increasing function defined in Ref. 



0. 



(5) 



Note that, in our 



analytical model, we will use this concept in order to introduce a minimum conductivity 
in the graphene sample, where the function $ p will be modeled by a free parameter to fit 



the experimental data and will take into account not only the carriers generated by the 
impurities but also other kind of phenomena that could contribute to induce carrier density 
The shear viscosity rj, in Eq. ([2]), can be calculated using 

where C v ~ 0(1) is a numerical coefficient, a = e 2 /ehc is the effective fine structure constant, 
e the relative dielectric constant of the substrate, and M the number of species of free 
massless Dirac particles js, 7|. Additionally, the entropy densities can be calculated according 
to the Gibbs-Duhem relation e + p = Ts. These equations have been derived under the 
assumption \u\ < c, and therefore the relativistic correction term, oc dp/dt, can be neglected, 
so that the classical Navier-Stokes equations are recovered. Note that, despite the high speed 
of the electrons, \u\ ~ 0.1c, the Reynolds number remains moderate, due to nano-metric size 
of the samples and the high kinematic viscosity of the electronic fluid in graphene. 



Kinematic Viscosity 



Based on Ref. [3j, the dynamic viscosity of graphene in a sample of linear size L , is given 
by Eq. ([6]). This equation can be written in the following form: 

M (k B T\ 2 h 

where we have introduced the characteristic frequency ujf = c/L , and by solving the appro- 
priate quantum Boltzmann equation, it is concluded that C v ~ 0.449. Eq. (J7J) can also be 
rewritten as rj = C v qJ 2 h/ Lq, where qf = huj /(k B T). Note that, in order for a classical (non 
quantum) picture of electron fluid to apply, the energy of excitations must be much lower 
than the thermal energy, i.e. 9^ < 1, the so-called collision-dominated regime. Taking a 
typical set of parameters (in MKS units), c = 10 6 , L = 10~ 6 , T = 100K, and 77/s ~ 0.2H/ks, 
we obtain rj ~ 10 -20 . Since the Reynolds number is dictated by the kinematic viscosity of 
the fluid, v, rather than by the dynamic one, 77 = pis, with no need of involving the mass 
density, it is of interest to estimate the kinematic viscosity of the electron fluid in graphene. 

To this purpose, we appeal to the definition of the Reynolds number as given in Ref. 
namely: 

s/k B k B Tu L 

Re = ' ( 8 ) 

rj/h nujf uq 
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where vq = cLq. By writing Re = uqLq/v and equating with the above, we obtain 

hu f s/k B 
k B T rj/s 

Using rj/s = js| and qj ~ 0.07, we obtain v ~ 10~ 2 . To be noted that, in spite 

of its extremely low dynamic viscosity, the kinematic viscosity of graphene is about four 
orders of magnitude larger than that of water. These four orders of magnitude are more 
than compensated by the large speed of the electrons, which is ultimately responsible for 
the sizeable values of the Reynolds numbers which can be achieved in graphene samples at 
micron scales. For instance, by taking u = 0.1c ~ 10 5 m/s, for a sample of 1 micron in 
length, we obtain Re ~ 20 for the global sample, and about Re ~ 0.04 on the scale of the 
impurities. 

Analytical Model Description 

In this work, we will treat impurities as circular rigid obstacles of diameter d. This choice 



is not arbitrary, but it is based on the fact that some experiments [ill, [15|, [16| suggest the 
the main scattering mechanism in graphene is due to strong neutral defects, with a range 
shorter than the Fermi wavelength, inducing resonant scattering. Thus, the diameter d can 
be interpreted as the characteristic length for the range of the interaction electron-impurity. 

Let us now assume that the electronic fluid moves in the x direction as a consequence 
of an applied electric field E, and V • (p r u) — (incompressible limit). Therefore, we can 
calculate the force F acting on a single impurity due to the electronic flow, as F = § II ■ dl, 

•o- 

where II is the stress tensor defined by 11^ = p5ij + p r UiUj + mj, with the viscous tensor 
7Tjj = t](dui/dxj + duj/dxi), and I is a unit vector perpendicular to the impurity circle. Here 
the indices i,j = 1,2 denote the coordinates x and y. Thus, solving the equations for the 
Dirac electron fluid, in the steady state (all time derivatives are neglected), and following 



an analogous procedure as in classical hydrodynamics 2l|, we obtain for the drag force, 
F d = E x , 

F D = X V v , (10) 

where A is a dimensionless parameter that depends on the Reynolds number. Here, v is the 
velocity of the fluid very far from the impurity. We first simulated single impurities with 
different sizes and different fluid velocities (see Fig.[TJ), obtaining that a linear approximation 



is justified in the range of relevance to this work. Note that, for a perfect fluid (v = 0), 
there is no drag force. However, from the point of view of the kinematic viscosity, electrons 
in graphene are far from being a perfect fluid. Therefore, we expect them to play a crucial 
effect on the drag force controlling the total conductivity of the sample. 

Let us denote by the impurity fraction, namely the ratio between the area not occupied 
by the impurities and the total area of the system, <fi = 1 — Nird 2 / (4A), with N the number 
of impurities in a sample of area A. We can thus propose the relation between v and the 
total current density as J = <ppv. Note that (pp is the existent electronic charge density in 
the graphene sample, since the volume fraction 1 — <fi is occupied by the impurities. 

Let us consider a representative elementary area of the sample of length d x in the direction 
of the flow and d y across it. With M circular impurities in this area, we can write M = 
4(1 — (f))d x d y / (ltd 2 ). The total force acting on the electronic fluid due to impurities (equal 
to the total force acting on the impurities due to the fluid), is given by F tot = MF^,. Here, 
the distance between impurities is taken sufficiently large enough to prevent the flow close 
to an impurity from affecting the flow nearby another impurity. 

In order to describe correctly the physics of graphene, we need to include in our model 
an extra feature. Due to the linear Dirac-Weyl spectrum of graphene, and the non-existence 
of a gap between the conduction and valence bands, the slightest amount of impurities 



or an external potential will induce charge carriers in the graphene sample [22], |23|, see 
Eq. (jSj). Thus, in our model, the total amount of carriers induced by the impurities will 
be proportional to the impurity concentration, (1 — (f))A, leading to an extra carrier density 
in the fraction of the sample occupied by the electronic fluid, (pA.The extra carriers are 
then given by 7(1 — 4>)/4>, where 7 is the proportionality constant that characterizes the 
impurity-electron interaction. 

Summing the forces, namely the Lorentz and drag forces, acting on the elementary area 
leads to: (ppEd x d y + 7(1 — <p)Ed x d y /(f> — F tot = 0, and by inserting the value of F tot , we 
obtain 

pud 2 ( (j) 7 



D 



71 E ■ (11) 



i 

Replacing Eq. ( fTUI) . taking into account that J = <ppv and Ohm's law, we can identify the 
conductivity as: 



a v \ Vi-0 > ' ^r^ + ^n , (12) 



p 2 ird 2 ( 4> 2 7 

a 



where we have introduced the coefficients <To = p 2 nd 2 / (4r/A), and a m i n = (?ol/p- This 
equation represents the key result of our paper ._An analogous derivation, for fluid dynamics 
in disordered media, can be found in Ref. 



24 



25| . Note that cr also can be written 



as o"o = n(e 2 ird 2 c 2 /12/c#Tz/A), where n = p/e is the electronic number density and v the 
kinematic viscosity. From this expression, we can see that the conductivity in graphene 
depends linearly on the carrier density, thus implying a constant mobility p = a/ne, in 
agreement with experimental observations jl[ Q]. In addition, our model can also explain 
why the mobility remains almost constant in the range of temperatures where v oc T" 1 
(see Ref. (9]), and the presence of a minimum conductivity in graphene, second term on 
the rhs of Eq. ( I12p . a m i n = eird 2 c 2, y /12k bTuX, which is independent of the carrier and 
impurity densities. Indeed, this model cannot explain, as other theoretical models, the 
sublinear behavior of the conductivity for the zero range impurity because, in that range, 
point defects and boundary conditions start to be dominant. In addition, it cannot describe 
the electron-phonon interaction either, since these have been excluded at the outset. All 
features above will make the object of future extensions of this work. 



Dimensionless Numbers 

For the numerical validation, and in order to obtain general results, we will use dimen- 
sionless numbers. For this purpose, we can rewrite Eq. ([3]) alternatively as 

« + (S. V ) 3+ 1 VP + -?* - JLVHJ = —E, (13) 

Ol p r p r C z Ot p r p r 



and therefore we obtain 



— + (u- V)u+-Vp+— 1 £-vV 2 u= , (14) 

at p r p T c l at p r 

where v represents the kinematic viscosity. Let us define the following relations: u = uqu', 
t = t^, (x,y,z) = L (x', y', z'), p r = p r oPU P = PoP'i an d E = E E', where the prime 
quantities are dimensionless, and Uq, t , L , p r0 , p , and E are characteristic values for 
the respective physical quantities. Thus, using the state equation e = 2p, we can deduce 
P — \Proulp' and e = |p r0 MQe'. Replacing these relations in Eq. ( Fl4|) . multiplying by and 
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using the relation u = we obtain, 

% + • V) u' + Ivy 

1 Tig ^ V 2 . Pq-EqLq p' A 

H -7; V w = 5 -C/ . 

3 c 2 at 1 u L p r0 UQ p' r 

To simplify this equation, we can equal the characteristic velocity to the Fermi speed, u = c. 
Therefore, we obtain 

§ + (u'-vv + -Lvy 

Pr (16) 
. d P' v w/2 PqEqLq p -., 

+ Ta77 r~ v n = 2 — r ■ 

We can identify two characteristic dimensionless numbers. The first one is the Reynolds 
number, which is, Re = ^^j^, and the second one, which we call 11 £ number" is defined by 
g _ p E Q Lp _ go_o ^gj-g y — _E L is the characteristic electrical potential. Using these 

PrO u Q PrO u Q 



definitions, we get 



9u' _,. 1 



+ l?f_ ±v , 2it = E l e , (17 » 

3 <9t' i?e p' r 
Note that this equation is dimensionless and therefore the universal features of the dynamics 
of the system are controlled only by the numbers Re and S: the latter measures the strength 
of the electric drive, while the former scales inversely with the dissipation opposing this drive. 
For notational simplicity, we will remove primes, leading to 

§ + ( 3 ■ v) u + J- v P + ff - i-v 2 a = e > J, (is) 

ot 3p r 3 at Re p r 



Numerical Results 



Fig. [2] illustrates the speed of the fluid for two different impurity densities, dark and 
yellow colors denoting low and high speeds respectively. An electric field of 1.77 V/m was 
applied in x direction (from left to right). Here we can see that for high impurity fraction 
(see Fig. |2} top), the speed of the fluid presents fluctuations on larger scales affecting the 
total conductivity of the sample. From the calculation of the electric current density and 
the electric field, we obtain the Ohm's law, giving a slope a. 

10 



The conductivity a is calculated from the numerical slopes and plotted as a function of 
the impurity fraction. The inset of Fig. [3] reports the comparison between the analytical 
solution, using Eqs. ( II 2p . showing an excellent agreement with the numerical data. For 
the fitting parameters, we obtain cq = (9.9 ±0.1) x 10~ 2 e 2 /h, and u m in — 3.4 ± 0.6. Note 
that there is a difference between our analytical model and the numerical simulations for the 
minimum conductivity. This difference is due to the fact that, for high impurity densities, the 
flow around one impurity starts to affect the flow around the others, and therefore, Eq. (]12p 
needs some additional terms. In particular, the approximation F tot = NFd does not hold 
anymore and non-linear correction terms should be included. Thus, while the minimum 
conductivity given by the analytical model is ~ 3Ae 2 /h, the simulations give ~ 4e 2 /h. We 
have verified that the conductivity of graphene, as computed in our model, does not show 
any appreciable dependence on the size of the system. 

In order to compare with experiments, we express the conductivity in terms of the ratio 
n/rii, where rti = (2.91 x 10 16 m _2 )(l — (f)), in our case. According to this expression and 
setting n = no, we obtain that £ = (l/a)n/rii = (1 — with a = 4.85 x 10~ 3 . Inserting 
this result into Eq. (|T2|) . we obtain a/ao = £(1 — l/£) 2 - Note that for values £ ^> 1, i.e. 
n/rii^> a = 4.85 x 10 -3 , this equation tends to 

a ~ <7 £ + Crmra = ~ V O min ■ (19) 

a Ui 

This corresponds to the linear dependence obtained by different theoretical models for 



graphene [13| 



flfl 



23 



22H3Q] 



In Fig. [31 we see the dependence between the conductivity and the ratio n/rij, and we 
clearly observe the prediction for the minimum conductivity of our model. The experimental 



3JJ , and compared with the results of the present 



data have been taken from Refs. 



work, showing good agreement. In Fig. [31 we also compare with the model proposed by 



Hwang et al. 



131 ] . where the impurities are located in a plane (substrate) parallel to the 



layer of graphene, with a separation 5 between the layers. In Fig. HI we compare our results 



with Coulomb impurity charges in random phase approximation (RPA) 13j. In the RPA 
model, the Boltzmann transport equation is used with impurities that are located randomly 
in the graphene sample. Our model shows good agreement in the slope with the RPA model, 
however, we achieve higher values due to the shift made by the minimum conductivity. 

11 



DISCUSSION 



We have developed an analytical model which accounts for a linear behavior of the 
conductivity with the electron density n, as well as with the ratio n/rii, in the limit 
n/rii^> 4.5 x 10~ 3 . In addition, it can also model the minimum conductivity in graphene as a 
consequence of the carrier density induced by the presence of impurities. Our model is based 
on a hydrodynamic description of electron flow in graphene, whereby Coulomb interactions 
are included through the viscosity of the electron fluid, and is valid in the collision-dominated 
regime. In this model, the impurities are treated as hard- sphere obstacles submerged on the 



electronic fluid, based on the fact that some experiments [111 Il5l . Il6| suggest that strong 
short-range neutral scatterers are the main scattering mechanism in graphene. Although this 
idea and the one about the long-range Coulomb scatterers are still object of controversy, the 
fact that the present analytical model can account for the conductivity of graphene suggests 
that indeed the short-range scattering models might be appropriate for graphene. 



This wor 
in Ref. 



£ is based on the hydrodynamic description of electrons in graphene proposed 



7|, which is a model developed for undoped graphene that neglects the electron- 
impurity and electron-phonon interactions. Here we have -extended- this approach by adding 
the electron-impurity interactions through a macroscopic porous media approach. Since this 
approach rests on basic conservation laws, it is supposedly very robust and independent 
on the validity of an underlying quantum Boltzmann equation, so long the microscopic 
interactions justify the build-up of a macroscopic viscosity (no superconductivity or other 
macroscopic quantum effects of that sort). Thus, our model is able to reproduce experimental 
results to a satisfactory degree of accuracy. 



For the set of parameters investigated in the present work, linear Ohm's law appears 
to apply throughout. However, based on Ref. Jq], by increasing the size of the impurities 
(less screening), non-Ohmic behavior could occur, due the onset of pre-turbulent phenomena 
within the graphene sample. It would be very interesting to verify such possibility by future 
experiments, as well as the inclusion of the electron-phonon interaction to model both, 
suspended and supported samples, at higher temperatures. 
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METHODS 



For the simulation, we use the hydrokinetic fluid solver proposed by Mendoza et al. 
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341 ]. adapted to two-dimensional flow in graphene 9[. The simulation was implemented on 
a grid of size 256 x 512 cells, representing a rectangular graphene sample of size 1.5 x 3/iin. 
We set up samples with a fixed number of impurities located randomly on the grid, each 
impurity covering one cell size, varying between 0.4 and 0.999. The Dirac-quasiparticle 
fluid in graphene has a kinematic viscosity v = 8.57 x 10~ 3 m 2 /s (see Ref. 9J), and by 
taking the Fermi speed no = 10 6 m/s as a characteristic speed, we obtain a Reynolds number 
Re = 350. Equating Re for both systems, in physical and numerical units, the cell size and 
time step are fixed to Sx = 5.86nm and 5t = 5.86fs. For a given temperature, T = 100K in 
our case, we can calculate the carrier density induced by the thermal energy with Eq. (Hj), 
no = 1.411 x 10 14 m -2 and therefore, using the approximate relation e = 2rinkBTn jsl, Q], the 
energy density e = 3.90 x 10 _7 J/m 2 and the density p r0 = 5.84 x 10 _19 kg/m 2 . Using the 
electron charge, we obtain the charge density, po — 2.26 x 10~ 5 C/m 2 . In numerical units, 
these values correspond to no = 4.845 x 10~ 3 , e = 2/3, and p r o = Po = 1-0, where the charge 
of the electron is e = 2.064 x 10 2 . Using the characteristic velocity Uq, we can calculate 
the value of the characteristic current density J = po«o = 22.6A/m or J = pou = 1.0 in 
physical and numerical units, respectively. On the other hand, to obtain realistic values of S, 
we use a characteristic electric field of E = 4.41V/m, which in numerical units corresponds 
to En = 10~ 6 . In this work, £ takes values from 10 -5 to 2 x 10 -4 . 



To model the extra carrier density induced by the impurities, as described in Eq. ([5]), we 
introduce an extra density charge Ap localized on each impurity position. Therefore, each 
impurity contributes a quantity p* = Ap dx 2 /A to the total charge of the sample, such that 
P = Pth + P*N = pth + Ap rii, where rii denotes the impurity density. This linear dependence 
between p and n^ is in qualitative agreement with experimental data 23J. We made several 
simulations for different values of Ap, finding that Ap = 60 leads to a minimum conductivity 
of 4e 2 /h. In our analytical model, this value corresponds to 7 = 60, in numerical units. The 
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simulations ran up to 5 x 10 5 time steps. 
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Figure 1. Drag force acting on a single impurity. Drag force Fq acting on a single impurity 
as a function of the graphene flow drift velocity for different impurity diameters. The solid lines 
represent the linear dependency of the drag force with the velocity of the fluid. In the inset, the 
dependence of the dimensionless parameter A on the impurity diameter is shown. 




Figure 2. Speed of the electronic flow. Absolute value of the velocity in graphene with 
multiple impurities, for two different impurity fractions, 0.952 (bottom) and 0.999 (top). The 
electric field is applied in the x direction (from left to right) and set up to 1.77 V/m. 
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Figure 3. Comparison between our analytical model and experimental data. Comparison 
between our results (stars) and experimental data for the conductivity a, as a function of n/rij. 
Data from Ref. [l| 



and from Ref. 



are represented by up and down triangles, from Ref. 



23J by circles and squares, 



311 ] by diamonds and crosses, for electrons and ho 



es respectively. Solid lines from 

n 

bottom to top, theory for separations 5 = according to Ref. [13[, our results, and theory for 
5 = 0.2nm according to the previous reference. In the inset (top), we show the conductivity as a 
function of 2 /(l — <fr), with an inset to observe the minimum conductivity. In the inset (bottom) 
we amplify the region close to the Dirac point. 
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Figure 4. Comparison between analytical approaches. Conductivity a as a function of n/rii 



for different types of scattering models 



131 ] . RPA is the conductivity calculated by using a random 



phase approximation with Coulomb scatterers. The unscreened Coulomb interaction would yield a 
conductivity smaller than the minimum value in graphene, over the entire range of gate voltages. 
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